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ABSTRACT 

We present the results of a study with the Swift Burst Alert Telescope in the 14 - 195 
keV range of the long-term variability of 5 low mass X-ray binaries with reported or 
suspected super-orbital periods — 4U 1636-536, 4U 1820-303, 4U 1916-053, Cyg X-2 
and Sco X-l. No significant persistent periodic modulation was detected around the 
previously reported periods in the 4U 1916-053, Cyg X-2 or Sco X-l light curves. The 
~170 d period of 4U 1820-303 was detected up to 24 keV, consistent with variable 
accretion due to the previously proposed triple system model. The ^46 d period in 4U 
1636-536 was detected up to 100 keV, with the modulation in the low and high energy 
bands found to be phase shifted by ~180° with respect to each other. This phase shift, 
when taken together with the near-coincident onset of the ~46 d modulation and the 
low/hard X-ray state, leads us to speculate that the modulation could herald transient 
jet formation. 
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1 INTRODUCTION 

Long-term "super-orbital" periodic variability has been ob- 
served in soft X-rays from over 30 X-ray binaries with the 
All Sky Monitor (ASM) on the Ross i X-ray Timing E xplorer 
(RXTE) satellite (see for example ISood et all [20071 '). Also 
known as 'long' or 'third' periods, these modulations are 
defined simply as any periodic variability greater than the 
orbital period. This very general definition of super-orbital 
variability is consistent with a number of physical mecha- 
nisms for the variation. Unlike binary orbital or neutron star 
spin periods, the mechanisms behind super-orbital variabil- 
ity are not well understood. The precession of a warped ac- 
cretion disc modulating the X-ray flux from the compact 
object is currently the favoured model (see for example 
IWiiers fc Pringld Il999l ; IClarkson etHI |2003al h This warp 
can be cause d by the gravitationa l field of the donor star (see 
for example lLarwood et al.lll996l h or by radiation pressure 
exe rted by the central X-r a y emitting regions (see fo r exam- 
ple IwHelSZlEi^ showed 
that while radiation driven warping gives a coherent picture 
of super-orbital periods, no t all systems agree with its pre- 
dictions (e.g. 2S 0114+650: iFarrell et al.ll200sh . 

Other models proposed to explain super-orbital vari- 
ability include variations in the accretion rate (caused by 
stellar pulsations of the donor star, irregular Roche lobe 
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overfl o w or as y mmet ric stellar wind capture; IWeiler et al.l 
1 19921 : iPsaltisj | 2006), precession of the neutro n star 
spin-axis in an X-ray p ulsar (|Trumper et al.l 1 19861 : 
IStairs. Lvne fc Shemarll200oh. variations in the location of 
an a ccreting hot-spot l|Rutten. van Paradiis fc Tinbergenl 
1992) and the presence of a third body in the sys- 
tem (i.e. a triple syste m; iPaul. Kitamotofc Makinol I2OO0I : 
IChou fc Grindlavll200ll h 

While most super-orbital periods are at best quasi- 
periodic, four X-ray binaries have been shown to exhibit 
stable super-orbital modulation over lon g timescales - Her 
X-l, LMC X-4, SS 433 and 2S 0114+650 (|Paul fc Kitamotd 
120021 : ISood et all l2006l h Other systems such as SMC X- 
1 have super-orbital periods which have be en observed to 
vary s moothly within a well-defined range (|Clarkson et al.l 
I20 038J ). while other super-orbital periods appear to evolve 
stochastically. It is thus clear that there is more than one 
mechanism behind the super-orbital modulations. 

A number of X-ray binaries were regularly monitored 
in hard X-rays by the Compton Gamma Ray Observa- 
tory (CGRO) Burst and Transient Source Experiment 
fBATSE) for the fir st ~1,600 days of the RXTE mission 
dHarmon et all 120041 1. Since the decommissioning of CGRO 
in June 2000, the lack of high energy response of the RXTE 
ASM has precluded long-term studies of super-orbital vari- 
ability simultaneously at low and high energies. The launch 
of Swift in 2004 saw the arrival of a wide field monitoring in- 
strument capable of complementing the RXTE ASM in hard 
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X-ray s - the Burst Alert Telescope (BAT; iBarthelmv et al] 
20051). In this paper we present the results of a search for 
the previously reported super-orbital modulation of 5 low 
mass X-ray binaries (LMXBs) in the BAT 14 - 195 keV 
light curves. 



2 SOURCES 
2.1 4U 1636-536 

4U 1636-536 was discovered following a series of type 
I X-ray bursts observed w ith the SAS 3 satellite 
(Hoff man. Lewin fc Doty 1977) and later classified as an 
atoll source ijHasinger fc van der Klidll989h . The orbital pe- 
riod of 3.8 h r was derived from optical photometr ic mea- 
surements by(Pedej^gji J _vaiiFjaradijs_&_ Lewin ( 1981). A dis- 



tance of 4.3 ± 0.6 kpc was derived by IJonker fc Nelemans! 
|2004h using type I photospheri c radius expansion bu rsts. 
However, more recent analyses bv lGallowav et all (|2006l ) de- 
rived a distance of 6.0 ± 0.5 kpc for a canonical neutron star 
(i.e. Mats = 1.4 M G R N S = 10 km). 

ICasares et al.l l|2006r i performed phase-resolved VLT 
spectroscopy of the irradiated donor in 4U 1636-536, using 
Doppler imaging of the NIB A4640 Bowen transition lines 
to determine a new set of spectroscopic ephemerides. A re- 
analysis of burst oscillations using these ephemerides allowed 
them to place the first dynamical constraints on the masses 
of the binary companions, constraining the mass function 
(f(M) = 0.76 ± 0.47 M ), mass-ratio (q ~ 0.21 - 0.34) and 
syst em inclinat i on an gle (i ~ 36 - 60°). 

IShih et al] (|2005l ) reported the onset of a large ampli- 
tude, statistically significant periodic modulation at ~46 d 
in the RXTE ASM data, interpreting this periodicity as due 
to an accretion rate variability related to the X-ray irradia- 
tion of the disc. These same authors noted a gradual decline 
in the X-ray flux of this source coincident with the appear- 
ance of the ~46 d variability, arguing that 4U 1636-53 was 
transiting from activity to quiescence d ue to a reduction in 
the mass accretion rate. iFiocchi et all (|2006l ). however, ob- 
served a hard tail dominating the spectrum above 30 keV 
using data from the INTEGRAL satellite, speculating that 
this high energy emission could be linked to the onset of 
relativistic outflows (i.e. jets). 

The hard state of atoll sources has been previously 
considered to be associated with the presen ce of a s teady 
jet and radio emission (see for example iFenderl 120061 ; 
Migli axTfc Fendeill2006l ). When compared w i th bla ck hole 
systems containing jets, iMigliari fc Fender! (|2006l ) found 
that the neutron star systems were less radio loud for a given 
X-ray luminosity, but that the jet power in both classes 
of object scales linearly with accretion rate. Neutron star 
LMXBs in the low/hard state are thus expected to contain 
weaker radio jets than their black hole counterparts. 

A number of radio observations of 4U 1636-536 were 
performed before the onset of the 46 d modulation, the de- 
creased X-ray flux, and the hardening of the X-ray spectrum, 
showing an absence of significant radio emission (see § 6.2.4). 
No radio observations targeting this source after the onset 
of the modulation or the transition to the low/hard state 
have been reported in the literature. As such, the presence 
of jets in this system is yet to be confirmed. 



2.2 4U 1820-303 

4U 1820-303 was discovered near the centre of the glob- 
ular cluster NGC 66 24 in 1974 by the Uhuru satellite 
l|Giacconi et alj Il974h . The subsequent discovery of in- 
tense bursts associated with thermonuclear flashes from this 
source confirmed it as a low mass X-ray binary containing a 
neutron star, m aking this the first id entified source of type 
I X-ray bursts (|Grindlav et all 1 19761 ). A distance of 7.6 ± 
0.4 kpc has been estimated using the bursts as a standard 
candle (|Kuulkers et al.ll2003h . 

The system consists of a white dwarf and 
a neutron star orbiting their common centre 

of mass with an orbit al period of 11.4 min 

dRappaport et all Il987l; IStella. Priedhorskv fc White! 



Il987f ). lHasinger fc van der Klid (Il989l) later classi fi ed 4U 
1820-303 as an atoll source, with IZhang et ail (|l998l ) 
observing the presenc e of kH z quasi-periodic oscillations 
(QPOs). iBloser et all (|2000bl ) observed that the QPO 
frequency was strongly correlated with the position on the 
colour-colour diagram. This result (in combination with the 
observed saturation of the QPO frequency at high accretion 
rates) lead them to conclude that they were observing the 
last stable orbit in the accretion d isc. 

IPriedhorskv fc Terrelll (|l984bl ) noted a periodic ~170 d 
variability in the X-ray flux, over which timescale the source 
moves from the soft/banana state to the hard/island state. 
This quasi-periodic ~170 d variability was later proposed 
to be a result of variable mass accretion due to tidally in- 
duced mod ulation of the orbital ecc entricity b y a third mas- 
sivc body (|Chou fc Grindlavll200ll) . Recently, iTarana et alj 
l|2007h reported the detection of a hard X-ray tail above 50 
keV from this source in the island state, concluding that 
this emission is likely due either to emission of non-thermal 
electrons or thermal emission from plasmas with a relatively 

hig h temperature. 

IZdziarski. Wen fc Gierliriskil (|2007al ) performed an 
analysis of the super-orbital modulation using RXTE ASM 
data. These authors undertook a study of hierarchical triple 
system models yielding the required quasi-periodic eccen- 
tricity oscillations, finding that the resulting theoretical light 
curves matched well with the observed ones. Zdziarski et al. 
i|2007bl ) reported a dependence on the X-ray spectral state 
of both the depth and phasing of the orbital X-ray modula- 
tion. This led them to favour a model in which the observed 
X-ray modulation is caused by scattering in hot gas around 
an accretion-stream bulge at the disc edge; as the mass ac- 
cretion rate changes with the ~170 d cycle, so the size and 
location of this bulge varies. The average flux calculated over 
the course of the super-orbital variability was found to be 
compatible with the model of accretion due to the angular 
momentum loss via emission of gravitational radiation. 



2.3 4U 1916-053 

4U 1916-053 was disc overed in 1974 by the Uhuru satellite 
l|Giacconi et al.l 1 19741 ) and is the most compact X-ray bi- 
nary to exhibit periodic intensity dips. Type I X-ray bursts 
each lasting ~10 s have been observed fr om this source with 
a recurrence time scale of several hours (|Becker et al . 1977; 
iLewin et al.lll977|). A distance of 8.8 ± 1.3 kpc has been 
estimated by Jonker fc Nelemans! l|2004l ) from photospheric 
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radius expansion during the bursts. 4U 1916-053 is classified 
as an atoll source and displays a similar path in the X-ray 
color-color diagram to 4 U 1820-303, althou gh with signifi- 
cantly lower luminosity (Blos er et al.ll2000al ). 

The X-ray dips have been ascribed to obscuration of the 
X-ray emitting region by a bulge in the out er accretion disk 
wher e the accretion stream meets the disc (|White fc Swank! 
1982). The dip recurrence period of ^3000 s is the shortest 
amon g the known dipping sources l|Swank. Taam fc White! 
1 19841 ). The optical counterpart was discovered as a V = 
21 mag optical star m odulating at a period of 3027.5 s by 
iGrindlav et alJ l|l988h . This optical period is ~1% longer 
than the X-ray dip period, raising the question as to which 
of the values indicated the orbital period. 

A 3.9 d beat period betw een the o p tical and X- 
ray periods was reported by IGrindlav! l|l992l ). man- 
ifesting as a modulation o f the X-ray dip shape 
ijChou. Grindlav fc Bloseill200lf ). Timing analyses of RXTE 
PCA data and quasi-simu l taneo us optical data performed by 
IChou. Grindlav fc Bloserl (|200ll ) lead them to conclude that 
the X-ray dip period was indeed the binary orbital period, 
with the slightly longer optical period produced through 
beating between the orbital period and the 3.9 d period, 
which they ascribed to the precession period of the accre- 
tion disc. They also set an upper limit of 2.06 x 10 -11 s 
s _1 on the \P\, arguing that the inferred stability and the 
partial clustering of X-ray bursts preceding the X-ray dips 
continues to suggest that the system may be a hierarchical 
triple. 

A regular long-term periodicity of 199 d (and a possi- 
ble period of 40.6 d) was reported i n the X-ray flux using 
data f rom the Vela 5B satellite by (jPriedhorskv fc Terrelll 
1984a!) , which the authors a scribed to disc precession . Later 
analysis of the same data bv lSmale fc Lochneij (Il992l ) found 
that the 199 d modulati on was not ve r y stro ng, with a con- 
fidence level of just 80%. |5 omer et al ] (|200lh analysed ~3.7 
yr of RXTE ASM data, finding no signifi cant peaks ab ove 
the red or white noise levels at any period. ISimonl (|2005l ) in- 
vestigated the long-term X-ray activity of 4U 1916-053 using 
7 yr of RXTE ASM data, finding that a cycle of about 230 - 
250 d occurred only in a limited time interval, materialising 
as vari ations of the recurrence time and amplitude of out- 
bursts. ISimonl (|2005h also concluded that the movement of 
4U 1916-053 along the atoll path in the X-ray color-color di- 
agram during the outbursts implied that the observed long- 
term variability was caused by variable accretion onto the 
neutron star, and not by disk precession (which would im- 
ply occultation and/or absorption effect s). A more recent 
analysis of 8.5 yr of RXTE ASM data by I Wen et all (|2006l ) 
failed to detect any long-term periodic variability. 



2.4 Cyg X-2 

Cyg X-2 is one of t he six Galactic LMXBs tha t are 



|2004| ) respectively. A binary orbital period of 9.8444 d 
was deduced from the optical behavio ur of the donor 
star (Ca sares. Charles fc Kuulkersl Il998l ). A high energy 
tail has been observed i n the X-ray spectrum du r ing a 
number of observations llDi Salvo. Robba fc Stella 20021 ; 



iasingcr & van der Khslll989l). and 



Bvram. Chubb fc Friedman! (1966) 



classified as Z-sources ( 
was first discovered by 
with a sounding ro cket experime nt. Type I X-ray bursts 
were reported by Smalc (1998), confirming a neutron 
star as the compact object. Distances of 7.2 ± 1.1 
kpc and 13 ± 2 kpc have been estim ated from opti- 
cal observations ijOrosz fc Kuulke"rslll999 | ) and from photo- 
spheric radius expansion in the bursts Uonker fc Nelemansl 



iPiraino. Santangelo fc Kaaretll2002l ; lLavagetto et alj|2006l ). 
although hard X-ray emissi on does not appear to be a con- 
stant feature of this source (|Lavagetto et al.ll2006l ). 

A 77 d periodicity in data from the Ve la 5B 
satellite was reported by Smale fc Lochnerl l|l992h . 
IWiinands. Kuulkers fc Smalel ( 19961 ) reported the de- 
tection of a 78 d period in RXTE ASM data from the first 
160 d of the mission, and showed that their result was 
supported by archival data from the Ariel 5 and Vela 5B 
missions. However, analyses of archived data from the all 
sky monitors on the RXTE, Gtnga, Ariel 5, and Vela 5B 
satelli tes and the s canni ng modulation collimator of HE AO 
1 by iPaul et all (2000) found no evidence for a stable 
modulation in the X-ray emission. Strong peaks at 40.4 d 
and 68.8 d in the RXTE ASM data, and at 53.7 d and 61.3 
d in the Ginga ASM data were observed by these authors, 
leading them to conclude that the super-orbital variability 
in this system (if real) is aperiodic. Dy namic analyses 
erform ed on 7 yr of RXTE ASM data by IClarkson et al.l 
2003bl ) revealed complex quasi-periodic behaviour, with a 
highly significant periodicity at ~59 d in a 300 d section of 
the light curve, with significant recurring features at ~40 d 
appearing sporadically. 



2.5 Sco X-l 

The Z-source Sco X-l was the first X-ray binary dis- 
covered and is the brightest p ersistent extra-solar X- 
ray source ijGiacconi et alJ [l962). Although the system 
properties are consistent with it containing a neutron 
star with a weak magn etic field, no type I burst s have 
been observed so far (Mir abel fc Rodriguesl [2003). Ra- 
dio emission from the vicinit y of Sco X-l was first 
reported by lAndrew fc Purtonl (|l968l ). with later high- 
resolution radio interferometry unveiling compact jets 
through the for mation of twin radio lobes with var i- 
able morphology ijFomalont. Geldzahler fc Bradshawll200ll ). 
A distance of 2.8 ± 0.3 kpc was determined from 
trigo nometric parallax (Bradsha w. Fomalont fc Geldzahlerl 
1999). The orbital period of ~18.9 hr was originally dis- 
cov ered in 85 yr of archived pho tometric optical data 
by iGottlieb. Wright fc Lillen (l200lh. and more recently 
observ ed in X-rays by IVanderlinde. Levine fc Rappaportl 
(|2003h - 

A high energy tail above 30 keV in the X-ray spectrum 
was o bserved in CGRO OSSE data bv lStrickman fc Barretl 
(120001) and later verified in RXTE and INTEGR AL data 
bv lD'Amico et all l|200ll ) and lDi Salvo et all (|2006h respec- 
tively. A periodic m odulation at 37 d was reported by 
iPeele fc White] { 19961 ) in 9 months of RXTE ASM data, al- 
though this modulation has yet to be confirmed. 



3 DATA REDUCTION 

While awaiting gamma-ray bursts the Swift BAT instrument 
accumulates data in survey mode. A given source is within 
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the field of view of the BAT ~10% of the time, with a typi- 
cal integration time per pointing of around 1000 s. Standard 
software was used to extract and clean detector plane im- 
ages in four energy bands for each observation (typically a 
few thousand seconds, often comprising several shorter seg- 
ments) and to calculate the arrays of weights corresponding 
to the visibility of each source within the field of view by 
each detector pixel during that observation. A program de- 
veloped for the purpose by the authors was then used to fit 
simultaneously the intensities of the sources and a series of 
background model components (equivalent functionality to 
that provided by this program is now available in the latest 
public software, though the public software was not avail- 
able at the time the analyses were performed). Corrections 
were applied for earth occultation and the known energy- 
dependent variation of the instrument sensitivity with off- 
axis angle. 



4 DATA ANALYSIS 

The BAT 14 - 24 keV (A), 24 - 50 keV (B), 50 - 100 keV 
(C) and 100 - 195 keV (D) light curves were used for these 
analyses, rebinned at 1 d. The average sampling rate for the 
reduced light curves prior to rebinning varied from ~0.8 - 
12.5 d" 1 . For comparison, the RXTE ASM 1.5 - 12 keV 1 d 
average light curves for the same time span were analysed to 
cover the low energy range of the spectrum. The 1 d binning 
was chosen for convenience as we were interested only in the 
long-term variability of the target sources. Table [T] lists the 
number of data points and average sampling frequencies for 
the BAT and RXTE ASM light curves for each source. 

All data points were weighted using 
the modified weighting schem e developed by 
ICorbet. Markwardt fc Tuellerl l|2007l ). In this method, 
individual data points are given weights based on both the 
non-uniform measurement uncertainties and the intrinsic 
source variability. The resulting weighted count rate R'i is 
given by: 



R'i 



Ri 



[(M) 2 + Vs. 



(i) 



where Ri is the raw un-weighted rate, / is a correction factor 
to account for underestimation of the error on each point, 
and Oi is the uncertainty for the point. V s is the estimated 
variance due to source variability, given by: 



Vs = 



£f =1 (^-*i) 2 Ef =1 (M 



TV- 1 



N 



(2) 



where R is the mean count rate and TV is the number of data 
points in the light curve. 

The / factor was derived for each of the energy bands in 
the BAT and RXTE ASM light curves by fitting a constant 
to the light curves of a source expected to be stable. The 
/ factor was determined by taking the square root of the 
reduced \ 2 value obtained from fitting data from the Crab 
(R. H. D. Corbet 2007, private communication). For cases in 
which the V s parameter was negative (applicable to sources 
which are very faint in the given energy band), V s was set 
to zero and the meth od reverts to t he "simple" weighting 
technique proposed bv lScargld l|l989l ). 

Power spectra were generated for each of the weighted 



light curves using the fast periodogr am fasper subrou- 
tine of the Lomb-Scargle per iodogram (|Lombll 19751 : IScargld 
1 19821 : iPress fc Rvbickil ll989l V The over sampling and hi- 
frequenc y factors - which set the period range and res- 
olution (IPress fc Rvbickil fl989l) - were both set at 2 for 
each analysis run, allowing us to search for periodic vari- 
ability in the ~1.3 - 1000 d range. The Lomb-Scargle tech- 
nique allows us to sample below the average Nyquist pe- 
riod (with reduced sensitivity) due to the unevenly sam- 
pled nature of the data, without creating problems due to 
aliasing as would be seen in the analysis of an evenly sam- 
pled data set using standard fast Fourier techniques. The 
frequency range and binning were chosen to give good cov- 
erage of the range covered by known super-orbital periods 
while at the same time providing sufficient resolution in 
the power spectra. The 99% white noise significance levels 
were e stimated using Monte Carlo sim ulations (see for ex- 
ample |KongIChaner^Kuuikeri[r998 : ). The 99% red noise 
significance levels were estimated using the redfit sub- 
routing which fits a first-order autoregressive process to 
the time series in order t o estimate the red-noise spectrum 
jSchulz fc Mudelsedl2002h . 



5 RESULTS 
5.1 4U 1636-536 

Figure E]] shows the RXTE ASM and Swift BAT light 
curves for 4U 1636-536 while Figure IB II shows the power 
spectra. The previously reported 46 d periodicity is clearly 
detected in the RXTE ASM 1.5 - 12 keV data, and is also 
significantly present in the BAT 14 - 24 keV, 24 - 50 keV 
and 50 - 100 keV light curves. No significant modulation 
above 100 keV is detected in the 1.3 - 1000 d range. In or- 
der to investigate the likely origin of the 46 d modulation, 
the raw (i.e. un-weighted) RXTE ASM, BAT A, BAT B and 
BAT C light curves were each folded over the 46 d period, 
setting phase zero arbitrarily to the start of the BAT light 
curves at MJD 53360. The profile of the modulation was ob- 
served to shift in phase between the 1.5 - 12 keV and 14 - 
24 keV energy bands, so that the profiles above 14 keV are 
almost 180° out of phase with the lower energy modulation. 

It can be clearly seen in Figure lXTI that the period of the 
modulation is not stable, with the time between successive 
peaks varying over time. In order to test how the drift in the 
period affects the width of the peak in the power spectra, we 
simulated light curves with the same sampling as the RXTE 
ASM, BAT A, BAT B and BAT C light curves with a stable 
sinusoidal modulation with a period of 46 d. Power spectra 
were then generated and compared with those produced us- 
ing the real data shown in Figure [Bll In all cases the peaks 
in the real power spectra were marginally wider than the 
simulated spectra. While this discrepancy could possibly be 
linked to the change in the period, the difference could also 
be due to the effects of noise or statistical scatter. 

Analysis of the BAT light curves prior to rebinning 
into 1 d average data sets found an additional peak 
at ~77 above the 99% white noise significance levels. 



ftp: / /ftp. ncdc.noaa.gov/pub/data/palco/softlib/redfit/redfit_preprint.pdf 
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Table 1. Swift BAT and RXTE ASM light curve properties for each of the four sources analysed. 



Source 


RXTE ASM 


BAT A 


BAT B 


BAT C 


BAT D 
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f b 

t-samp^ 


N 


fsamp 


N 


fsamp 


N 


fsamp 


N 


fsamp 


4U 1636-536 


481 


0.77 


483 


0.78 


485 


0.78 


485 


0.78 


494 


0.80 


4U 1820-303 


452 


0.72 


419 


0.75 


405 


0.72 


426 


0.76 


429 


0.76 


4U 1916-053 


485 


0.78 


379 


0.69 


384 


0.70 


386 


0.70 


390 


0.71 


Cyg X-2 


570 


0.91 


504 


0.81 


493 


0.79 


507 


0.81 


507 


0.81 


Sco X-l 


485 


0.77 


483 


0.77 


470 


0.75 


466 


0.74 


465 


0.74 



a number of points in the light curve 

" average light curve sampling frequency in units of d 
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DC 
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Figure 1. The light curves of 4U 1636-536 folded over the 46 d 
period. From bottom: RXTE ASM 1.5 - 12 keV profile, BAT A 
14 - 24 keV profile, BAT B 24 - 50 keV profile, and BAT C 50 
— 100 keV profile. The vertical dashed lines indicate the phases 
corresponding to the minimum and maximum of the RXTE ASM 
profile. The count rates are normalised to the mean of each light 
curve and the scale for each profile is the same. 



iFarrell, O'Neill fc Soodl l|2005t ) have shown that by impos- 
ing a known periodic sampling rate, any spurious peaks in 
the power spectrum produced by aliasing or spectral leakage 
should be shifted to other frequencies. By simulating light 
curves with periodic modulations at 46 and 77 d with the 
same sampling as the BAT data we found that rebinning 
the light curves should increase the significance of a real 
peak. The significance of the 77 d peak was dramatically 
reduced in power spectra of the rebinned BAT light curves 
(see Figure IB1[I and is well below the dominant 99% white 



noise significance levels. We thus identify the 77 d peak as 
spurious. 



5.2 4U 1820-303 



Figure [£J shows the RXTE ASM and Swift BAT light 
curves for 4U 1820-303 while Figure IB2I shows the power 
spectra. The previously reported ~170 d periodicity is 
clearly detected in the RXTE ASM 1.5 - 12 keV data as 
a single broad peak, and is also significantly present in the 
BAT 14 - 24 keV light curves. No peaks around 46 d are 
significantly present in the 24 - 195 keV power spectra. The 
1.5 - 12 keV and 14 - 24 keV light curves were thus epoch 
folded over the 170 d period, setting phase zero arbitrarily to 
the start of the ASM light curve at MJD 53353. The count 
rates were normalised to the mean rate (Figure [2]). Unlike 
4U 1636-536, no phase shift was observed between the low 
and high energy profiles. 

The period of the super-orbital modula tion of 4U 1820- 
303 h as been known to vary over time (Zdziarski et al. 
l2007al ). As such, to test the effect of a varying period on 
the width of the peak in the power spectra, we simulated 
light curves with a stable sinusoidal modulation at 170 d 
with the same sampling as the RXTE ASM and BAT A 
light curves. We then generated power spectra for the simu- 
lated light curves and compared them with the power spec- 
tra shown in Figure IB2I As with 4U 1636-536, the peaks 
in the real power spectra were marginally wider than those 
generated using the simulated data. While this discrepancy 
could again be linked to a changing period, there is no more 
than a suggestion of possible slight broadening of the peaks. 

The simulations clearly show that the width of the peak 
is a function of both the period of the modulation and the 
number of cycles in the data set. Therefore, in order to com- 
pare the width of the super-orbital peak in the 4U 1820-303 
power spectra with that of 4U 1636-536 this relationship 
needs to be taken into account. We thus took the full width 
of the peaks at 50% of the maximum power, normalised 
it by the period of the modulation and then multiplied by 
the number of cycles in the light curves. Comparing the nor- 
malised peak widths shown in Table[2] we see that the values 
for the real peaks are comparable for both 4U 1636-536 and 
4U 1820-303, and are consistently marginally greater than 
the width of the simulated peaks. 

In addition to the ~170 d modulation, a peak at around 
~100 d marginally exceeds the 99% white noise significance 



6 S. A. Farrell et al. 



BAT B 




0.5 1 1.5 2 

PHASE 

Figure 2. The light curves of 4U 1820-303 folded over the 170 d 
period. From bottom: RXTE ASM 1.5 - 12 keV profile, BAT A 
14 - 24 keV profile and BAT B 24 - 50 keV profile. The vertical 
dashed lines indicate the phases corresponding to the minimum 
and maximum of the RXTE ASM profile. The count rates are 
normalised to the mean of each light curve and the scale for each 
profile is the same. 



Table 2. Comparison between the normalised widths of the 
super-orbital peaks in real and simulated power spectra for 4U 
1636-536 and 4U 1820-303. 



Data Set 


4U 1636-536 


4U 1820-303 




real 




real 




RXTE ASM 


0.94 


0.82 


l.n 


0.90 


BAT A 


1.07 


0.88 


0.94 


0.86 


BAT B 


0.97 


0.88 






BAT C 


1.00 


0.88 







level in the 24 - 50 keV power spectrum. However, the lack 
of a significant peak at this period in the other BAT power 
spectra (in particular between 14 - 24 keV) leads us to con- 
clude that this peak is spurious. No other significant modu- 
lation above 50 keV is detected in the 1.3 - 1000 d range. 

5.3 4U 1916-053 

Figure \Mi shows the RXTE ASM and Swift BAT light 
curves for 4U 1916-053 while Figure IB3I shows the power 
spectra. No significant modulation in any of the light 
curves is apparent in the 1.3 - 1000 d range, indicating 
that either the 199 d and 40.6 d periodicities detected by 




-1.5 1 » .'-U-IJ-M-I 

-3 -2.5 -2 -1.5 -1 -0.5 

LOG(FREQUENCY) [LOG(D 1 )] 

Figure 3. Log-log power spectra of the 1.5 - 12 keV RXTE ASM 
light curves of 4U 1916-053 covering the lHomer et all J200lh time 
frame (grey solid lines) and the time frame covered by the BAT 
data (black solid lines). The red noise models for each power 
spectrum are plotted as dashed lines for comparison. 

IPriedhorskv fc Terrelll (|l984al ) were spurious, the mecha- 
nism for these modulations was not active during the time 
of our observations, or that the signal was not sufficiently 
stable and/or persistent over the entire observing window 
for it to be detected in a single power spectrum. 

The power spectrum of the RXTE ASM data for 4U 
1916-053 appears to be dominated by white noise. This re - 
sult is in contrast to the results of Iflomer et al.l (|200 ll ). 
who found that red noise dominated the power spectrum 
of the first four yea rs of RXTE ASM data. Analysis of the 
iHomer et all l|200ll ) data using our analysis methods con- 
firm that red noise was indeed dominant in the RXTE ASM 
data, but has since dropped in significance. Figure[3]presents 
a comparison between the powe r spectra gen e rated for the 
RXTE ASM data covering the IHomer et ail (|200ll ) epoch 
(MJD 50150 - 51500) and the time span used for our anal- 
yses (MJD 53360 - 53900), clearly showing the drop in red 
noise. 



5.4 Cyg X-2 

Figure EJ] shows the RXTE ASM and Swift BAT light 
curves for Cyg X-2 while Figure [B4l shows the power spec- 
tra. A strong, broad peak at ~190 d is present in the RXTE 
ASM 1.5 - 12 keV and BAT 14 - 24 keV power spectra. 
These peaks are located at approximately 0.5 yr, a com- 
mon location for spurious peaks (see for see for example 
iFarrell et al.l [20051 ) . No other significant peaks are present 
in any of the energy bands. A peak at ~40 d is present in 
the power spectrum from the RXTE ASM data , similar to 
the 40.4 d signal detected by iPaul et ail l|2000l ). Although 
the power of this peak exceeds the 99% white noise signifi- 
cance level, it is well below the red noise significance level. 



5.5 Sco X-l 

Figure EH shows the RXTE ASM and Swift BAT light 
curves for Sco X-l while Figure IB5I shows the power spec- 
tra. No periodic variability greater than the dominant 99% 
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white noise significance level is present in the RXTE ASM 
or BAT power spectra in the 1.3 - 1000 d range. 



6 DISCUSSION 

The typical shape of the X-ray spectrum for LMXBs con- 
sists of a thermal Comptonisation component (sometimes 
empirically modelled as a power law with exponential cut- 
off) curtailed at low energies by photoelectric absorption, 
often with the addition of a low-energy black body com- 
ponent arising either from thermal emission from the inner 
regions of the disc or very close to the surface of the ne utron 
star (see for example iBloser et al.ll2000bl ; lBarredl200ll. for a 
review). In the low/hard state an additional power law com- 
ponent is sometimes req uired in order to mo del the hard X- 
ray tail (see for example iTarana et al.| [2007). with the hard 
X-ray emission attributed variously to Comptonisation by 
high-temperature thermal or non-thermal electrons in the 
accretion disc corona or in compact jets. 



6.1 The ~170 d Period in 4U 1820-303 

The ~170 d period in 4U 1820-303 has been attributed to 
variable accretion linked to the presence of a third body 
in this system, whereby the eccentricity of the inner or- 
bit is modulated by tidal effects produced by th e outer 
body (|Chou fc Grindlavll200ll ; IZdziarski et al.ll2007al ). Tran- 
sitions between the high/soft and low/hard states over this 
timescale are linked to variability in the accretion rate 
llBloser et al . 2000b) and the movement in and out of the 
inner region of the accretion disc (see for example iMalzad 
120071 ). thus modulating the thermal Comptonisation region 
of the spectrum. This thermal Comptonisation component 
has been seen to dominate the spectrum up to ~50 keV, 
above which 4U 1820-303 has only been detected on one 
occasion during the low/hard (island) state, during which 
a high energy non-thermal tail was found to dominate 



|Tarana et allbOOTf ). 

The ~170 d period has been seen to correspond to nor- 
mal atoll source motion in the colour-colour diagram, with 
the minimum of the mo dulation coinciding with the transi- 
tion to the island state (IBloser et al. 2000b). If the appear- 
ance of the high energy non-thermal tail is regular and tied 
to the island state (and thus the minimum of the super- 
orbital modulation), a hardening of the spectrum would be 
anticipated every ~170 d. We would therefore expect to de- 
tect the modulation above 50 keV, anti-correlated with the 
low energy modulation. If the hard tail in 4U 1820-303 pref- 
erentially appears on a regular basis during the transition 
to the island state, we would also expect the low and high 
energy modulations to be out of phase. However, as men- 
tioned in § 5.2, no such phase shift is apparent in the folded 
profiles (Figure [2}. 

The super-orbital period should be present at least up 
to the 50 keV threshold, as a varying accretion rate would 
be expected to modulate the entire thermal Compton region 
of the spectrum. To test whether 4U 1820-303 would be de- 
tected by the BAT above 24 keV, we simulated spectra (1000 
s integration times) using the BAT response files for the 
various island state spectral models applied to RXTE and 



INTE GRAL data by IBloser et~ai] i|2000bh and lTarana et all 
|2007j). respectively. 

Flux upper limits during the minimum of the ~170 
d modulation (i.e. the island state) in the three BAT en- 
ergy bands a bove 2 4 ke V were calcula t ed fo r each of the 
IBloser et all |2000b) and ITarana et all (|2007l ') island state 
spectral models. Table [3] lists the derived flux upper lim- 
its in the BAT B, C and D energy bands fo r each of the 
IBloser et ail l|2000bl ) and ITarana et all (|2007h models. The 
projected flux in each band for each spectral model was also 
calculated, with the normalisation parameters set so that 
the 14 - 24 keV model count rate equalled the measured 
rate in this band (the only BAT band in which we can be 
certain 4U 1820-303 was detec ted). 

In all cases except the Tarana et al.l l|2007l ) simple 
Comptonisation and the IBloser et al.l (l2000bl ) cut-off power 
law + disc black body models, the expected flux was far in 
excess of the derived upper limits in the 24 - 50 keV and 
50 - 100 keV bands. Thus, within the uncertainties of this 
method, the non-detection of the ~170 d modulation above 
24 keV is mo st consistent with th e sof test island state spec - 
tral models in lTarana etaH l|2007l ) and lBloser et all (|2000bh . 



6.2 The ~46 d Period in 4U 1636-536 

6.2.1 Formation of Compact Jets? 

iFiocchi et al.l (|2006l ) linked the recent reduction in X-ray 
flux and the appearance of a high energy tail in the spectrum 
of 4U 1636-536 to the onset of jet formation. Interestingly, 
the 46 d modulation first appears in the RXTE ASM data 
around the same time as this transition to the low/hard 
state. Figure [4] shows the coarse binned RXTE ASM light 
curve up to MJD 53900 for 4U 1636-536 with the signal- 
to- noise (S/N) of the 46 d peak (as taken from overlapping 
power spectra generated from 600 d light curve sections with 
a 100 d overlap). The significance of the modulation can 
clearly be seen to increase markedly with the decrease in flux 
(the significance of the 46 d periodicity did not exceed the 
99% significance levels until a S/N of ~7 was reached). The 
appearance of a high energy power law tail in the spectrum 
in conjunction with the first detection of the 46 d period 
and the reduction in flux implies a link between the th ree 
phenomena. If the assertion made bv lFiocchT et all (|2006l ) is 
correct (that the transition to the low/hard island state is 
linked to jet formation), the 46 d modulation could be tied 
to transient jet formation, or possibly to the precession of 
jets (linked in turn to the precess ion of a warpe d disc) in a 
similar fashion as seen in SS 433 (|Margonlll984T l. 

As mentioned in § 5.1, the ~46 d modulation is clearly 
seen in each of the folded profiles shown in Figure [1] with 
the phase of the modulation shifting with increasing energy 
so that the BAT C (50 - 100 keV) profile is almost 180° out 
of phase with the RXTE ASM (1.5 - 12 keV) profile. A sim- 
ilar phase shift was observed between th e RXT E ASM and 
INTEGRAL light curves by IShih et al.l (|2005l '>. which was 
attributed to the high and low energy modulations origi- 
nating at different locations. In comparison, the profiles of 
the light curves of 4U 1820-303 when folded over the 170 d 
period show no such phase shift (Figure [2}. 
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Table 3. Derived flux upper lim its for 4U 1 820-303 in th e island state i n the BAT B, C and D energy bands 
using the model parameters in lBloser et al ] |2000bh and lTarana et al ] j2007lV 



BAT Energy Band 




Flux Upper Limits (x 10 11 erf 


;s cm 2 s 1 ) 






CPL + BB» 


CompTT + BB^ CPL + DBB^ 


CompTT£ 


CompTT + PL« 


B (24 - 50 keV) 


5.5 


5.5 5.2 


4.8 


5.5 


C (50 - 100 keV) 


1.5 


1.5 1.6 


1.6 


1.7 


D (100 - 195 keV) 


7.4 


7.0 7.6 


7.9 


8.6 



cut-off power law plus black body model jBloser et al.ll2000bl) 
thermal Comptonisation plus black body mod el llBloser et al. 2000 
cut-off power law plus disk black body model IjBloser et al.ll2000 

thermal Comptonisation model (Tarana et al. I l2007l) 

thermal Comptonisation plus power law model ifTarana et al.l 2007) 




I 

49500 50500 51500 52500 53500 

Time (MJD) 



Figure 4. Average RXTE ASM 1.5 - 12.0 keV count rates 
(crosses) from 4U 1636-536 compared with the signal-to-noise of 
the 46 d peak in the power spectra (solid line). The arrows indi- 
cate the times of two previous radio observations made with the 
ATCA, while the dashed line indicates the timing of the INTE- 
GRAL observations which first detected the high energy tail in 
the X-ray spectrum. 

6.2.2 Interpreting the Phase Shift 

The detection of the 46 d modulation in both the high and 
low energy realms of the spectrum indicate that both the 
thermal and non-thermal components are varying over this 
timescale. The phase shift of the modulation is inconsis- 
tent with the standard warped precessing disc model for 
super-orbital variability, as periodic obscuration or reflec- 
tion would be expected to cause the high and low energy 
modulations to vary in phase. We interpret this phase shift 
as an indication that the 46 d modulation arises from two 
separate (yet closely linked) phenomena. If the low energy 
modulation seen in the RXTE ASM light curves is indeed 
due to variable obscuration of the central emission region 
(either the inner regions of the disc or the neutron star it- 
self), and the high energy modulation seen in the BAT data 
is a result of a varying viewing geometry of the beamed 
emission from the jet, it is possible to see two similar mod- 
ulations in different energy bands that are out of phase. In 
this scenario, the jet column will be pointing towards us at 
roughly the same time as the bulge in the disc is obscuring 
the central emission region. Thus we would see an increase 
in the high energy flux at the same time as a decrease in the 
low energy flux (as observed in the RXTE ASM and BAT 



data), with the high and low energy flux varying smoothly 
as the disc and jet precess together. 

Many systems which show super-orbital modulation as- 
cribed to disc precession do not show any evidence for jet 
formation, and in some cases jet formation is inhibited due 
to the presence of a highly magnetised neutron star (e.g. Her 
X-l). However, if a warped precessing disc is present, any 
jets which form would expect to be coupled to the warped 
disc and should themselves precess. The jet would be aligned 
with the normal axis of the inner disc, so the scenario in 
which the bulge is phased to obscure the inner disc at the 
precise time as the jet points towards us is unlikely. 



6.2.3 Transient Jet Formation? 

A somewhat simpler explanation would be that the 46 d 
modulation is linked to periodic transitions between the 
high/soft and low/hard states, leading to transient jet for- 
mation. In this alternate scenario, the inner region of the 
disc moves in as the accretion rate intensifies, resulting in 
an increase in the flux of seed photons and therefore an en- 
hancement in the thermal Comptonisation (i.e. soft) compo- 
nent of the spectrum. In this high/soft state the conditions 
are not believed to be sufficient for jet formation: the disc 
is thought to be geometrically thin yet optically thick, with 
a relatively small corona as most of the accretion power is 
believed to b e dissipated in the disc itself (see for example 
lMalzadl2007l '). The hard region of the spectrum during this 
state is characterised by a weak, steep power law, generally 
interpreted as inverse Compton up-scattering of soft pho- 
ton s by non-the rmal electrons in the corona (see for exam- 
ple Erret] |200j) . As the accretion rate decreases, the inner 
regions of the disc move further out and the disc cools, be- 
coming geometrically thick yet optically thin. In this state 
most of the accretion power is believed to be transported 
into the corona, which in turn is thought to power the com- 
pact jet. The soft thermal component from the disc decreases 
as the flux of seed photons lessens, and the strong hard non- 
thermal component appears in conjunction with the jets. 
Thus an ant i- correlation between the hard and soft X-ray 
flux would result, with both fluxes modulated over the ~46 
d period. In this scenario the 46 d period would represent 
the characteristic timescale between disruption and recovery 
of the disc. 
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6.2.4 Searching for Radio Emission 

Recently, iMigliari et all (|2007f ) showed there to be a direct 
link between radio jet emission and the X-ray state in the 
neutron star X-ray binary GX 17+2. They also showed that 
the presence of a hard X-ray tail was coupled to the low/hard 
X-ray state, implying a link between the presence of radio 
jets and the appearance of the hard tail in the spectrum. 
Thus, if transient jet formation linked to spectral state tran- 
sitions is the cause of the 46 d period, we would expect radio 
emission from this source to be modulated in phase with the 
hard X-ray modulation, as the jets form an d radio emission 
ensues (see for example ICorbel et al.l[200^ ). 

A review of the Australian Telescope Compact Array 
(ATCA) online archive (ATOA) shows that two previous 
radio observations of 4U 1636-536 were performed on 1995 
September 2 and 1998 April 4, with upper limits of 0.12 mjy 
and 0.15 mjy derived a t 4.8 GHz and 8.6 GHz respectively 
l|Berendsen et al.ll2000h . Both observations were performed 
during the high/soft state, when compact jets and resulting 
radio emission would not be expected. More recent obser- 
vations were performed with the ATCA on 2007 June 29, 
30 and July 1 as part of a multi-wavelength collaboration 
to study LMXBs of all kinds at various epochfl While this 
data is still proprietary, the principal investigator of these 
observations was kind enough to communicate the prelimi- 
nary results; a 3<r upper limit of ~0.6 mjy was derived at 
8.6 GHz (J. Dickey 2007, private communication). 

The 2-10 keV X-ray flux in the low/hard state mea- 
sured with INTEGRAL was ~1 x 10~ 9 erg cm" 2 s* 1 
(M. Fiocchi 2007, private communication), giving a lumi- 
nosity of ^2 x 10 36 erg s _1 for a distance of 4.3 kpc 



I Jonker fc Nelemansl l2004h and ~4 x 10 36 erg s -1 for a 
distance of 6 kpc l|Gallowav et al"1l2006l ). According to the 
X-ray/radio luminosity correlation for neutron star LMXBs 
derived by Migliari & Fender (2006), 4U 1636-536 should 
thus have a radio luminosity (from measurements around 
8.5 GHz) of somewhere between (3 - 20) x 10 28 erg s" 1 . 
The 8.5 GHz flux density would thus be ~1 mjy, barely 
sufficient for a detection to be made in the recent ATCA 
observations. 

Comparing the timing of these recent observations with 
the RXTE ASM 1.5 - 12 keV and the BAT 14 - 50 keV 
light curves (Figure [SJ shows that they occurred just after 
the maximum of the modulation in the soft X-ray band, and 
just before the increase in hard X-ray flux. These observa- 
tions thus appear to have been performed before the tran- 
sition to the low/hard state and the predicted phase of jet 
formation. Further support for this theory is obtained when 
the timing of the INTEGRA L observations — du ring which 
the hard tail was detected bv lFiocchi et alj l|2006j ) — is com- 
pared with the soft X-ray RXTE ASM light curve (Figure 
It can clearly be seen that the detection of the hard tail 
occurred during the minimum of the soft X-ray modulation, 
at the same time as we would have expected to see a peak in 
the BAT hard X-ray light curve had monitoring data been 
available during this epoch. 
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Figure 5. RXTE ASM (top) and BAT 14 - 50 keV (bottom) light 
curves covering MJD 54220 - 54420, showing the timing of the 
recent ATCA observations of 4U 1636-536 (arrows), which took 
place just prior to the transition to the low/hard island state. The 
dashed horizontal lines indicate the average BAT count rate and 
the average high/soft state count rate in the RXTE ASM. The 
latter line shows that the source re-enters the high/soft state at 
the peak of the ~46 d modulation in the soft X-ray band. 
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Figure 6. The RXTE ASM light curve of 4U 1636-536 cover- 
ing MJD 52650 - 52750, showing the timing of the INTEGRAL 
observations (arrow) which detected the hard X-ray tail. These 
observations took place during the drop in low energy flux, which 
presumably coincided with a rise in the high energy flux. As the 
observations took place after the decommissioning of CGRO and 
prior to the launch of Swift, long-term monitoring data in hard 
X-rays is not available during this epoch. 



6.3 Non-detections 

No evidence for the previously reported super-orbital peri- 
ods in 4U 1916-053, Cyg X-2 or Sco X-l was found in the 
RXTE ASM or any of the BAT light curves. The power spec- 
tra were instead dominated by low-frequency noise, with no 
peaks that would indicate the presence of real long-term pe- 
riodic modulation. The modulations for 4U 1916-053 and 
Cyg X-2 were or iginally reported using Vela 5B data, with 
later analyses by I Wen et all (|2006l ) of 8.5 yr of RXTE ASM 
data finding no significant stable modulation from either 
source. The detection of super-orbital modulation in Sco 
X-l was reported using only the first 9 months of RXTE 
ASM data. It is thus apparent that either the reported de- 
tections were spurious, the mechanisms behind the super- 
orbital variability in each of these three systems has since 
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ceased, or the modulation does not persist over a sufficient 
fraction of the observing window to be detectable by our 
techniques. Unfortunately, the BAT light curves do not cover 
a sufficiently long epoch for us to subdivide the light curves, 
thus precluding dynamic analyses to search for signals that 
are potentially transient or variable in frequency. 



7 CONCLUSIONS 

We have analysed the RXTE ASM and Swift BAT light 
curves in the 1.5 - 195 keV range for 4U 1636-536, 4U 
1820-303, 4U 1916-053, Cyg X-2 and Sco X-l. We found no 
evidence for periodic variability at the previously reported 
super-orbital periods from 4U 1916-053, Cyg X-2 or Sco X- 
1 in either the RXTE ASM or BAT data. In contrast, the 
~170 d period in 4U 1820-303 is clearly detected up to 24 
keV, consistent with the modulation arising from variable 
accretion due to tidal effects from the third body in the 
triple system. The 46 d period in 4U 1636-536 was signifi- 
cantly detected in the RXTE ASM data and up to 100 keV 
in the BAT data, with the profile of the high energy mod- 
ulation found to be ~180° out of phase with that seen in 
the 1.5 - 12 keV band. This behaviour is inconsistent with 
the modulation arising solely from variable photoelectric ab- 
sorption due to the precession of a warp in the accretion 
disc, leading us to interpret it as evidence that the high and 
low energy modulations arise from two separate yet closely 
linked phenomena. The appearance of the 46 d period in 
conjunction with the transition to the low/hard state is in- 
triguing, leading us to speculate that this modulation could 
be tied to the formation of transient compact jets. In this 
scenario, the low energy modulation would result from a 
thermal Compton component of the spectrum which varies 
with the accretion rate. The high energy modulation would 
thus be tied to the appearance of the non-thermal hard tail 
arising from the formation of compact jets during episodes 
of decreased accretion. Further monitoring observations at 
radio wavelengths are being sought in order to confirm this 
conclusion. 
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Figure Al. 4U 1636-536: the 1.5 - 12 keV fiXTE ASM light curve (bottom panel) and the 14 - 24 keV, 24 - 50 keV, 50 - 100 keV and 
100 - 195 keV Swift BAT light curves (top panels). 
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Figure A2. 4U 1820-303: the 1.5 - 12 keV RXTE ASM light curve (bottom panel) and the 14 - 24 keV, 24 - 50 keV, 50 - 100 keV and 
100 - 195 keV Swift BAT light curves (top panels). 
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Figure A3. 4U 1916-053: the 1.5 - 12 keV RXTE ASM light curve (bottom panel) and the 14 - 24 keV, 24 - 50 keV, 50 - 100 keV and 
100 - 195 keV Swift BAT light curves (top panels). 
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Figure A4. Cyg X-2: the 1.5 - 12 keV RXTE ASM light curve (bottom panel) and the 14 - 24 keV, 24 - 50 keV, 50 - 100 keV and 
100 - 195 keV Swift BAT light curves (top panels). 
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Figure A5. Sco X-l: the 1.5 - 12 keV RXTE ASM light curve (bottom panel) and the 14 - 24 keV, 24 - 50 keV, 50 - 100 keV and 100 
- 195 keV Swift BAT light curves (top panels). 
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Figure Bl. 4U 1636-536: normalised power spectra of the 1.5 - 12 keV RXTE ASM light curve (bottom panel) and the 14 - 24 keV, 
24 - 50 keV, 50 - 100 keV and 100 - 195 keV Swift BAT light curves (top panels). The red horizontal dashed lines and blue solid 
lines indicate the 99% white and red noise significance levels respectively. The arrow indicates the location of the previously reported 
super-orbital period of 46 d. The power spectra were each normalised to the maximum power present, the 99% white or red noise levels, 
whichever was greatest. 




Figure B2. 4U 1820-303: normalised power spectra of the 1.5 - 12 keV RXTE ASM light curve (bottom panel) and the 14 - 24 keV, 
24 - 50 keV, 50 - 100 keV and 100 - 195 keV Swift BAT light curves (top panels). The arrow indicates the location of the previously 
reported super-orbital period of ~170 d. 
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Figure B3. 4U 1916-053: normalised power spectra of the 1.5 - 12 keV RXTE ASM light curve (bottom panel) and the 14 - 24 keV, 
24 - 50 keV, 50 - 100 keV and 100 - 195 keV Swift BAT light curves (top panels). The arrow indicates the location of the previously 
reported super-orbital period of 199 d. 
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Figure B4. Cyg X-2: normalised power spectra of the 1.5 - 12 keV RXTE ASM light curve (bottom panel) and the 14 - 24 keV, 24 - 
50 keV, 50 — 100 keV and 100 - 195 keV Swift BAT light curves (top panels). The arrow indicates the location of the previously reported 
super-orbital period of 78 d. 
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Figure B5. Sco X-l: normalised power spectra of the 1.5 - 12 keV RXTE ASM light curve (bottom panel) and the 14 - 24 keV, 24 - 50 
keV, 50 - 100 keV and 100 - 195 keV Swift BAT light curves (top panels). The arrow indicates the location of the previously reported 
super-orbital periods of 37 d. 



